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By Jack Funk and Richard H. Rhyne 


SUMMARY 


Vertical -tail loads were measured in turbulent air on a four-engine 
jet bomber. The results showed large and regtilar load oscillations which 
were lightly damped. Conqparlson of experimental results with discrete- 
load calculations indicated that discrete-gust calculations underestimated 
the loads by 30 to ^K) percent and gave no indication of the oscillatory 
characteristics or low danping. Calculations based on p,ower spectral 
analysis, on the other hand, reflected the general frequency character- 
istics of the measured loads and gave a better estimate of the tail loads. 
The present results strongly indicate that discrete-gust calculations 
for gust loads on vertical tails may seriously underestimate the gust 
loads for airplanes having lightly damped lateral oscillations. 


INTRODUCTION 


Reference 1 suggests that the gust loads on the vertical -tall sur- 
face can be calculated satisfactorily by the sharp-edge-gust equation 
with -no alleviation due to unsteeidy-lift effects or airplane motions. 

This recommendation was based on some early and limited flight tests of 
an XB-15 and an 0-2H airplane in the late 1930's. Since these early 
tests, changes in airplane configuration and the large increase in speed 
have served to conplicate the airplane responses to rough air. It was 
desirable, therefore, to reassess the problem of gust loads on the 
vertical -tail surface in order to determine whether the sinplified rela- 
tion of the sharp-edge-gust formula is still applicable. 

A fll^t investigation was undertaken with a jet bomber in order to 
obtain experimental data on the magnitude and characteristics of the load 
on the vertical tail in rough air and to assess the significant factors 
affecting the tail loads. Flight tests were made at two center-of- 
gravlty positions in order to provide a measure of the effects of changes 
in stability. The results of these test measurements are presented 
herein. Calculations of the loads on the vertical-tail surface due to 
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gusts were made as a supplement to the flight-load measurements in order 
to determine how accirrately they could be predicted by existing proce- 
diires. Two types of calculations were considered: One was based on 

discrete-gust techniques, as siiggested by reference 1; and the other, 
involving the statistical techniques of power spectral analysis, was 
based on continuous gust histories. The results of the calculations 
based on both methods are presented and compared with the flight-test 
measurements . 


SYMBOLS 


a lift-curve slope for wing, per radian 

a^ lateral acceleration, g un i ts 

a_ nomal acceleration, g units 


b wing span, ft 

Cq yawing-moment coefficient, N/qSb 



acn 


lateral-force coefficient, Y/qS 

lateral -force coefficient on vertical tail above strain 
meas\aring station 



“p ap 

Cy 


/vt 
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(%) 

c 

Db 

d 


S(Cy) 


vt 


g 


K, 


•g 

% 

Kz 


kz 

kg 

L 

^vt 


mean aerodynamic chord of wing, ft 

mean aerodynamic chord of vertical tail, ft 

differential operator with respect to 
damping coefficient 

acceleration due to gravity, ft/sec^ 
gust factor (ref. 3) 

nondlmensional radiun of gyration about principal longitudinal 
axis, 

nondimens lonal radius of gyration about principal vertical 
axis, kg/b 

nondlmensional radius of gyration about vertical stability 


. V 


p o P P 

axis, V K 2 cos t] + Kx sin'll! 


radl\is of gyration about principal longitudinal axis, ft 

radius of gyration about principal vertical axis, ft 

variation of lift coefficient when penetrating a sharp-edge 
gust, expressed as a fraction of final lift 

scale of t\rrbulence, ft 

tail length (from center of gravity of airplane to 25 percent 
mean aerodynamic chord of vertical tail), ft 


m 


mass of aiip)lane, W/g, slugs 



h 
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n 

q 

r 

S 


% 


T(oi) 

t 

U 


V 

Ve 


W 

Y 




P 


^^g 

P 

Po 


yawing moment on airplane, ft-lb 
dynamic pressure, Ib/sq ft 

yawing angular velocity, xadians/sec ' 

wing area, sq ft 

nondlmenslonal time parameter based on span, ^ 

frequency-response function 
time, sec 

true gust velocity, fps 

derived gust velocity, fps (ref. 3) 

true airspeed, fps 
equivalent airspeed, fps 

weight of aiiTplane, lb 
lateral force on aiipjlane, lb 

aerodynamic load, on 96 -squeire-foot area of vertical tall 
above strain measuring station, lb (also used as subscript) 

angle of slldeslip, radians 

angle between principal longitudinal axis of inertia and 
flight path, deg 

relative-density coefficient based on span, m/pSb 
airplane mass ratio, 2W/paScg 
density of air, slugs/cu ft 

standard density of air at sea level, 0.002378 slug/cu ft 
spectrum of longitudinal conponents of atmospheric turbulence 
spectrum of lateral conponents of atmospheric turbulence 



mCA TN 37^1 


5 


i|r angle of yaw, radians (also used as subscript) 

fl reduced frequency, <o/V, radlans/ft 

(0 frequency, r adlans/sec 

£^3 undamped natural frequency, radlans/sec 

Dot over quantity denotes first derivative with respect to time. 

METHOD AND TESTS 


The method consisted of measuring the Irput-gust history and the 
load on the vertical tail during fli^t of the test aiiplane In rough 
air. Tvro measurements of the turbulence input were obtained. Derived 
gust velocities were obtained by using the normal acceleration to provide 
gust velocities for use In the discrete-gust calculations, and the rapid 
fluctuations of the eiirspeed were recorded to provide an input for the 
power spectral calculations. The loads on the vertical-tall surface 
were measiured with calibrated electrical resistance strain gages mounted 
on the three spars of the vertical tail. 

Two flights of approximately ^ miles eadh were made, one at a 
center-of -gravity .position of percent mean aerodynamic chord and one 

at a center-of-gravity position of 30.6 percent mean aerodynamic chord. 

The runs were made in clear rough air at an altitude of approximately 
1,500 feet above the terrain and at an airspeed of approximately 390 knots. 
The pilot was Instructed to use as little control motion as was consistent 
with safe flight. 


AIRPLANE AND HJBTRUMENTATION 


The aiiplane used in this investigation was a Jet-powered medium 
bomber. A line drawing of the airplane is shown in figure 1. The char- 
acteristics of the aiiplane pertinent to the present analysis are sum- 
marized in table I. 

Vertical -tail loads were obtained from electrical wire resistance 
strain gages mounted on the three spars and from lateral-acceleration 
measurements at the vertical tail. The strain gages were installed 
approximately l4 inches above the horizontal tail as shown in figure 1 
and were calibrated in a manner similar to that of reference 2. The 
various gage responses were not combined electrically but were recorded 
individually and coiriblned numerically in the calibration equation. The 
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strain-gage outputs were recorded on an l8-channel osclUograph loavlng 
elements with frequency responses that were flat to 6o cps or above. 

An MCA alr-dajnped recording accelerometer was installed near* the root 
of the vertical-tail surface (as shown in fig. l) in order to obtain 
the inertia load. The accelerometer had a natural frequency of l4.7 cps 
and was daaped to about 0.7 of critical. 

In order to obtain data for evaluation of derived gust velocities, 
two accelerometers were installed near the front spar at approximately 
the nodal points of the wing fundamental bending mode (spanwlse location 
approximately 24-5 Inches from fuselage center line). The accelerometers 
were the strain-gage type having a natural frequency of .about 12 cps and 
were oil-damped to 0.7 critical. Locating the accelerometers at the 
nodes eliminated the effect of the fundamental wing bending mode on the 
acceleration measiorements. The two acceleration measurements were com- 
bined electrically to eliminate the effects of roll and were recorded 
on a galvanometer with a natural frequency of about 12 cps. 

In addition to the measurement of turbulence obtained from the ver- 
tical acceleration, turbulence measurements were obtained from the 
recorded airspeed fluctuations. Inasmuch as the airplane is longitudi- 
nally in sensitive to turbulence at the higher frequencies, the airspeed 
fluctuations at these frequencies are a direct measure of the tiirbulence. 
For the present investigation, it was estimated that this condition would 
apply for the i*ange of frequencies above l/h cps. The airspeed fluctua- 
tions were measured with a standard MCA airspeed-altitude recorder which 
had. a natural frequency of about 100 cps. The pitot-static head was 
located on a boom about 6 feet forward of the nose to minimize errors 
in static presstire. The pitot lines were made as short as possible, 
balanced, and danped with orifices to provide a frequency response which 
was shown by calibration to be flat to about 10 cps. 

Measurements were made of the fuselage flexibility by using two 
yaw-attitude recorders, one mounted at the center of gravity of the air- 
plane and the other mounted in the tail. Differences in the yaw-angle 
reading of the two recorders gave the change in angle of sideslip of the 
tail due to fuselage flexibility. 

In addition to the above data, measurements were obtained of the 
control-surface positions, rolling velocity, pitch attitude, and yaw 
angle . 


EVALUATION OF DATA AND BISULTS 


As an indication of the general characteristics of the records, 
short sections of some of the pertinent quantities measured axe shown 
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in figure 2 . The methods used in the evaluation of these records can he 
conveniently separated into two parts; one pertaining to the discrete- 
gust procedure and the other involving power spectra as will he 
described. 

As a starting point for the discrete-gust and spectral evaluations, 
500-second sections of the pertinent records were read at 0.1-second 
Intervals along the time history. From the 0 . 1 -second-interval readings 
of strain and lateral acceleration, the vertical-tail airload 'vreis deter- 
mined from the relation 

Vertical -tail airload = Structural shear load + Inertia load 

where the Inertia load is the product of the mass of the structure 
extending above the strain-gage station and its measured lateral accel- 
eration. The evaluation of the structural shear load was made for incre- 
mental loadings from the ti*im load which was taken as the mean load for 
each test run. The shear load was determined from the shear strain gages 
located on the three spars. Static calibration of the shear gages indi- 
cated that the shear readings were affected by the bending moment in the 
structure so that the total shear at the root is given by the relation 

Structural shear load = ASg^^ + BSg2 + -*^1 ®^2 ^^3 

where 6g]_, 8g2^ and 5 g^ are the shear-gage responses and 8]32.^ ^2^ 

and are the b ending-moment-gage responses. The evaluation of sev- 

eral of the larger loads indicated that the contribution of the bending 
moment woiiLd add less than 1.5 percent to the shear load; therefore, the 
effect of the bending moment was neglected in the shear evaluations. 


Flexibility 

Bending of the fuseleige under the load on the vertical-tail surface 
causes a small change in the angle of sideslip at the tall which acts to 
reduce the load on the vertical tail. This change in sideslip angle due 
to fuselage flexibility was evaluated by using the two yaw-attitude 
recorders installed at the center of gravity and at the tail of the air- 
plane. Inspection of the yaw-attitude record at the tall indicated that 
the change in sideslip due to the fundamental fuselage mode was so small 
that it was within the reading error of the records. Thus, it appears 
that the effect of fuselage flexibility, if any, will be that due to 
"static" bending under the tail load. 

As a preliminary check on the magnitude of the effects of the static 
fuselage bending on the tail loads, incremental changes in yaw attitude 
corresponding to the change in vertical -tail airloeid from maximum to 



8 


MCA 5PN 37^1 


minimum or miniimmi to maximum (such aa shown hy the solid circles in 
fig. 2) were evaluated from hoth the center-of -gravity emd tail yaw- 
attitude records. The difference between the incremental change at the 
center of gravity and at the tail was tahen as an incremental angle of 
sideslip at the tail due to fuseleige flexibility. The associated change 
in load due to fuselage flexibility was then calculated from the conven- 
tional lift equation hy use of the Incremental sideslip data. The results 
of this check Indicated that, on the average, fuselage flexibility reduced 
the loads on the vertical-tail surface by about 10 percent. 


Discrete-Gust Analysis 

Peak vertical -tail loads .- Q!he frequency distributions of the peak 
tail loads were evaluated from the time-history plot of the 0.1- second 
readings of vertical -tail airload as follows: The mean load was first 

determined from the 0.1-second readings. Then the values of the peak 
positive and negative load increments between successive crossings of 
the mean load line were tabulated into 100-pound class intervals. The 
distributions obtained for the t^ro center-of -gravity positions are given 
in table II. 


From the tabulated data, the average flight miles required to exceed 
given values of load increment M^AIvt) were determined by the fol- 
lowing relation: 



Total miles 

N(AIvt) 


where is the number of peak values exceeding a given value of 

load increment. The results obtained for the two test center-of -gravity 
positions are shown in figure 3 and represent the basic description of 
the load history. It should be noted that these data include the effects 
of turbulence variation between test runs. 


Derived gust velocity .- Distributions of derived gust velocity were 
calculated from the peak nodal accelerations according to the derived- 
gust -velocity formula from reference 3: 


®%iax 


apoSVe%e 

2W 




where 

®'%iax ®J^lane m a ximiTm nondimensional normal acceleration in g units 
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The nodal accelerations were evaluated from the steady-flight position 
of the acceleration trace for the largest peak "between consecutive 
crossings of the steady-flight trace position. Derived gust velocities 
were then calculated from these peak readings and are ta"bulated into fre- 
quency distributions in table IH. These distributions were used to 
determine the average flight miles required to exceed a given gust veloc- 
ity in a manner similar to that used for the vertical -tail loads, and the 
resiilts obtained are plotted in figure 


Power Spectral Analysis 

Vertical -tail load spectra .- In addition to the discrete-gust 
analysis, power spectra of the gust loads on the vertical-tall surface 
were obtained to show the frequency content of the load histories and 
to be compared with the theoretical calculations. The load, spectra 
were calculated by the procedure suggested by Tukey in reference 4 and 
reviewed in reference 5* Initially, calculations were made to cover 
the range from 0 to 5 cps eind were based on a 3.> 000-point san^le taken 
at 0.1- second intervals along the time history. Autocorrelation func- 
tions of 4o lags were made from which 40 power estimates were obtained 
over the frequency range of 0 to 5 cps. These results are shown in 
figures 5(a-) and 5(b) for the tests at center-of -gravity positions of 
24.3 and 30.6 percent mean aerodynamic chord, respectively, and are 
designated by the term "wide filter." (This numerical calculation can 
be viewed as equivalent to scanning the true spectrum with a roughly 
triangular filter having a l/2-cps base width. ) 

Because of the peakedness of the measiured spectra, a second anal- 
ysis was made to determine more precisely the shape in the vicinity of 
peak power. The second spectra were obtained from 600-point samples 
taken at l/2-second intervals along the time histories. (Only every 
fifth reading of the original data was used. ) Autocorrelation functions 
of 60 lags were calculated from which 60 power estimates covering the 
frequency range of 0 to 1 cps were obtained, (if 6o power estimates 
had been obtained for the frequency range of 0 to 1 cps for all 3^000 
readings at a 0.1-second interval, 300 lags for the autocorrelation 
function would have been required. The procedure used permitted a more 
practical calculation without any appreciable loss in reliability for 
the estimates over the desired range.) The effective filter width for 
this analysis is about l/l5 cps. The results are also shown in fig- 
Tires 5(Q’) and 5(t>) and are designated by the term "narrow filter." 

Because of the limited record lengths, individual power spectral 
estimates are not too reliable - roxighly, *30 percent of the values shown 
in figxire 5- The reliability of the overall power for each case is, of 
course, considerably greater, with the root-mean-square values estimated 
to be reliable to about ±10 percent of the value obtained. 
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Gust spectra.- 3ja order to ceilctilate tlie spectra of tlie vertical- 
tall load for comparison with the spectra of the measured load, it was 
necessary to obtain a measure of the spectra of the lateral gust veloc- 
ity. Although no direct measurements of the lateral-gust-velocity 
spectra were made, references 6 and 7 suggest that a satisfactory meas- 
\ire may he obtained from the power spectra of the longitud inal gust 
velocity on the assun^jtlon that atmospheric turbulence is isotropic. 

The longitudlnal-gust-velocity spectra were, therefore, evaluated from 
0.1-second-lnterval readings of the airspeed records. 

The power spectra of the airspeed fluctuations were, characterized by a 
very large peak in power at the low frequencies. The effects of such large 
peaks is to introduce some errors into the spectral analysis at the higher 
frequencies. In such cases, it is frequently helpful, as pointed out in 
reference 8, to reduce the power at the peaks by properly filtering the 
original data. This procediire is sometimes referred to as "prewhltenlng . " 

The power at the low frequencies i-ras reduced in aii 5 )litude before 
the spectral analysis was made by applying a high-pass filter to the 
data. The filter consisted of siibtractlng a moving average over a 
5-second Interval from the time history of airspeed. This technique is ' 
discussed in reference 8, in which it is shown that this procedture alters 
the power spectra by the factor 



where T is the total time interval over which the moving average is 
taken. For T of 5 seconds, almost all the power below 0.2 cps is 
removed. The spectra obtained from the filtered data were multiplied 
by the reciprocal of the foregoing filter in order to congpensate for 
the initial filtering. The resulting airspeed spectra are presented in 
figure 6. As a simple measure of the relative Intensity of the turbu- 
lence for the two test runs, the root-mean-square values of the spectra 
above 0.2 cps ane given in figiire 6. The power below 0.2 cps was not 
included because little power appears in the load spectra below 0.2 cps. 


CALCULATIONS 

Discrete-Gust Calculations 


The simplest approach for calcxilating the load on the vertical tail 
for a discrete gust is to con^iute the load resulting from a steady-state 
angle-of-sideslip change due to the gust. This procedure has been sug- 
gested for calculating gust load on a vertical -t^l surface in reference 1 
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and wltli some modification is the "basis for most ciirrent design proce- 
diires. The load on the vertical tail is obtained by siibstituting U^e/^ 

for the angle of sideslip due to the gusts so that 

’'vt (1) 

where ^vt is the load on the 96-square-foot area above the strain 

measuring station, and jCVo) is the slope of the lateral-force coef- 

\ P/vt 

flcient for the 96-square-foot area. By use of the preceding equation, 
the average flight miles to exceed given values of peak vertical -tail 
load were computed from the faired distributions of the derived gust 
velocity of figure i)-. The values obtained are presented in figure 5 as 
the curves labeled "discrete ceilculatlons . " These calculated loads are 
for the 96-square-foot area above the strain measuring station so that 
they can be compared directly with the measured loads. The value 
of ] used for the calculations was obtained from the fli^t- 

\ P/vt 

test results of reference 9* 


Spectral Calculations 

The basic relation between the power spectrum of the response of a 
linear system to a random disturbance and the spectrvan of the disturb- 
ance is given as (see ref. 5) 

$o(cd) = |T(cD)|%i(co) (2) 

where 

<C>o(^) power spectmmi of response or output 

T(o 3) frequency-response function 

<Ii^(a)) power spectnmi of disturbance or input 

and where bars designate the absolute value of the complex quantity. 

The application of this relation to the calculation of the vertical -tail 
load requires the determination of (l) the frequency-response function 
T(co) for the load on the vertical tall due to unit sinusoidal lateral 

gusts and ( 2 ) the lateral -gust -velocity spectnmi. The methods used in 

determining these two functions are considered in the succeeding section. 

Frequency-response function . - For the present analysis, the fol- 
lowlng assumptions were made: ("l) the airplane did not roll, (2) the 
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vertical-tail loads were due to side gusts only, and (5) the pilot had 
little effect on the loads. Inasmuch as the flight records indicated 
very small roll angles due to the action of the gusts, the restriction 
to the no-roll condition appeared Justified. The smal 1 roll angles 

^ahout li ^ also indicated that the asymmetrical component of the ver- 
tical gust velocity over the wing, which might give rise to sideslip motion 
and to load on the vertical tail, had little effect in the present case. 


Lag in lift was taken into consideration hy applying the gust pene- 
tration function k 2 (s) for a Mach number of 0.6 (from ref. 10) on the 

assuorption that the greater parts of the side-force anii yawing-moment 
coefficients come from the vertical-tail surface. This procedure is 
thought to give a satisfactory approximation to the effects of lag in 
lift. (Actually, the lag in lift ^s found to have very little effect 
on the calculated load spectra, the effect being a further attenviatlon 
of the spectra at the higher frequencies. However, the lag in lift did 
significantly affect the second moment of the spectra, which was needed 
in subsequent calculations of peak load. ) From the preceding considera- 
tions, the load on the vertical tail due to side gusts is given by the 
expression 

Y^(t) = qS^CYp^^p(t) 

The first term on the right-hand side of equation (5) is the load 
due to sideslip, the second term is the load due to yawing velocity, and 
the last term is the load due to the gust. The lag-in-lift function due 
to gust penetration kg is usually given in terms of airfoil chord 

length. It was changed to the time argument by the relation s = 

^vt 


+ qS 




iif(t) 


dU{ti) 

k2(t - ti)- ^ —- -dtp 


( 5 ) 


The frequency-response function for vertical -tall load 
obtained by substituting 



is 


U(ti) = U(cD)e^l 
•4r(t) = ilf(cD)e^^ 
P(t) = p(co)e^'t 

ixxst 


Yvt(-t) =Yvt(co)e' 
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into equation (5) and solving for - 777 — ?— • Uie result is 


u(cd) 






which may he written as 


03y^(£o) = qS^Cyr^^T^(cn) + ^ T^(ai) 


W 


[p(0,) + IQCo.)] 


i^) 


where 


j r»oo 

k2(t)e"^^dt 

0 

= © 

= IfS} 


The sideslip frequency-response function Tp(o)) and yaw frequency- 
response function T^(a)) may he determined from the conventional lateral 

equations of motion for two degrees of freedom such as are given in ref- 
erence 11. If the conventional stability axes shown in figure 7 
used, these equations become; 




Sideslip: 


Db - P + ^2pb Db - I ^ - 

cyr r% , . 


( 5 ) 
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Yaw: 


-CnpP - 





T]r = 



( 6 ) 


In equations (5) and ( 6 ), the operator Djj has been used to denote 

the differential — ^ where sj, is the nondimenslonal time parameter 

, ds^j “ b 

The lag-in-lift function which is usually expressed in terms of 

mean aerodynamic chord is converted to the new variable by the rela- 
tion b = 9*3cvt* two equations are solved simultaneously for T^(o>) 

and Tp(o3) . The results obtained are given by 


where 






P(<i>) + iQ((n) 
2 

V2 V 


Tp(03) 


CyqI) 


|p(cd) + iQ(cD^ 


-k£cjo? + A^im + B 
y2 V 


B = 


A - - ^p 

VbKzs^ 

*^p ^ ^r^p _ *^p^r 


( 7 ) 


(8) 
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Siibstltuting eg^tlons (7) and (8) into equation ( 4 a) gives the frequency- 
response function for the loads on the vertical tall as 


^vfcC“) 


where 



'^p(^r)vt , VCnr 

- (Oh + ^ i icD i03 

Vb^zs^ 

^(oi) + iQ(o))j 

Vl 

[cct)^ - + 2dicD] 

1 


( 9 ) 


2d 




2 ^ 

2t‘bKZ6^^ 


Estimation of stability derivatives .- The stability derivatives 
necessary for calculating the transfer functions were estimted in the 
following ma n ner; The coefficient for the 96 square feet of 

vertical -tail surface above the strain gages was obtained from the tests 
of reference 9. The coefficient Oy for the whole airplane was esti- 

mated by using the load measurements on the vertical tail and the lat- 
eral accelerometer at the center of gravity. The ratio of the total 
side load obtained from the lateral accelerometer to the measured load 
on the vertical tail was determined. Multiplying this ratio by (CV ^ 

\ ^pyvt 

gave Cy for the whole airplane. This method assumes that the effect 
P 

of flexibility on the fuselage aerodynamics is the same as for the verti- 
cal tail. The value of Cy obtained in this manner is thoiight to be 

P 

conpensated for the effects of steady-state flexibility. On the assTJnrp- 
tion that the yawing moment is produced by the forces on the verticsil 
tail, Cq was estimated from 
P 




where Oy ' refers to the derivative for the complete tail area 
P 

(127.8 square feet) and was obtained from (Gy \ on the basis that 

V Pjvt 

the derivatives were proportional to the ratio of the areas. 
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The values for Cy and 
■^r 



were estimated from the relation 


Qfr 



^vt 

b 72 


Cnr 



ivtf 

h2 




^vt 

vt b/2 


where the derivative 



refers to the 96 square feet of area 


above the strain-gage station. The stability derivatives estimated in 
the above manner are tabulated in table IV. 


By use of the stability derivatives in table IV and equation , 
the frequency-response functions for the vertical-tail load were calcu- 
lated and the arplitude squared of the frequency responses are shown in 
figure 8 for the two center-of -gravity positions. 


Input-gust spectra .- The spectra of the input-gust velocity were 
based on the spectra of airspeed fluctuations. However, as indicated 
prevloiisly, the airspeed fluctuations measure the longitudinal coniponent 
of the turbiolence, whereas the vertical-tail loads axe essentially the 
result of the side or lateral turbulence. On the assuiEption that atmos- 
pheric t-urbul6nce is isotropic, these two spectra are related. The spec- 
trum of the lateral con^onent of atmospheric turbialence, as indicated in 
reference 12, appears to be approximated very well by the following 
expression: 


* ^ 2 L 1 + 


( 10 ) 


where 


cfjj root-mean- square gust velocity 

The expression for the spectmm of the longitudinal component for iso- 
tropic txirbulence corresponding to equation (lO) is given by 


$x(«) 


C3 2L 1 


( 11 ) 


(See appendix of ref. 7.) If equations (lO) and (ll) are compared, it 
can be seen that the ratio of $y(ft) to depends upon the value 

of L HTiri varies with the value of fl. Reference 12, together with an 
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examination of the airspeed spectra of figure 6, suggests that L Is 
of the order of 1,000 feet or greater for the turbulence encotintered. 

For L of 1,000 feet and greater, $y(ft) is approximately 1-5 times the 
value of $x(^) values of greater than about 0.002, which in 

the present case corresponds to frequencies above 0.2 cps. The input- 
gust spectra for the calculation of the vertical -tail loads were accqrd- 
ingly obtained by multiplying the longitudinal spectra (fig. 6), obtained 
from the airspeed measurements, by 1.5- ■ 


The gust-load spectra for the vertical -tall surface were calculated 
by multiplying the input spectra discussed previously by the frequency- 
response functions (fig. 8) and are also shown in figure 5 for con^jari- 
son with the measured load spectra. 


Effect of errors in stability derivatives .- In order to evaluate 
the reliability of the calculation, the effects of errors in the stabil- 
ity derivatives on the calculations were considered. Uofortunately, 
there is no exact way to check the reliability of the stabilily deriva- 
tives . However, experience suggests that Gy and are reliable 

to ±10 percent and that 0^1^ is reliable to ±15 percent. The other 

stability derivatives have a much smaller effect on the results, and 
errors in these derivatives are of little importance. 


Calculations- were made to show the effect of errors in Cy , Cn ^ 

and Cnj. on the root -mean- square load. Each of these derivatives was 

individually reduced by 10 percent and a new value of the root-mean- 
square load was calculated by use of the gust-lpput spectra of figure 6 
and equations (2) and ( 9 ). The results showed '(that reducing both %: 

and by 10 percent reduced the root-m^an-square tall load about 


8 percent; a 10-percent reduction in Cn^ reduced the root-mean- sqTxare 
tall load less than 1 percent, and a 10-percent reduction in Cnp 


Increased the root-mean-square tail load by 4 percent. Therefore, the 
overall uncertainties in the stability derivatives are roughly estimated 
to yield uncertainties in the calculated root-mean-square loads of the 
order of ±10 to ±15 percent. 


Calculation of Peak-Load Distribution 

One of the principal objects of gust-load calculations is to esti- 
mate the magnitude and frequency of load peaks Tor use in design. If 
atmospheric turbtilence is considered to be a simple Gaussian disturbance, 
the average nuniber of peak loads per second exceeding a given value is. 
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as indicated in reference 12, given liy 



( 12 ) 


"where Oyyt root-mean-sguare tail load and $(cn) is the power 

spectrum of the tail load. 

A sinple evaluation of equation (12) for the calculated load spectra 
of figure 5 Indicated that the shape of the calculated curve did 

not conform to the measurements. The character of the difference was 
a less rapid increase in the calculated load "with increasing fli^t miles 
than that for the measured load. This characteristic difference be"tween 
flight measurements and calcialations based on Gaussian theory has been 
noted in other s"fcudies and is discussed in reference 12 in which it is 
attributed to 'a lack of homogeneity in the lntensi"ty of atmospheric tur- 
bulence. Reference 12 indicates that flight-test -data on peak gust loads 
can generally be synthesized by a combination of Gaussian turbiilence 
dis"turbances which "vaiy in intensi"fcy but which have the same spectral 
foim. For these conditions, the nuniber. of peak loads exceeding a given 
"vhlue of Yvt is given by 



where is the root-mean-square load for the individual Gaussian 

components. 

This model "was applied by determining the root-mean-square gust 
velocity for 50 segments of equal length for the present test data. The 
root -mean-square loads were then "taken proportional to the gust 

intensities for the corresponding segments. Equation (15) "was then 
evaluated for the data of both center-of -gravity positions and the 
res\ilts ob"tained are shown in figure 5* 
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DISCUSSION 

General ChEiracteristlcs of Measured Load 


[The general ctiaracteristics of tlie loads on the vertical tail in 
rough air are indicated by the sample load history shown in figure 2. 

The load history in rough air exhibits a continuous oscillation at a 
frequency of about 0.4 cps with varying amplitude. This load time his- 
tory is cheiracterlstlc of a lightly damped system responding to an irreg- 
ular but continuous exciting force. Comparison of the record samples 
shovm in figure 2 indicates that both the loads and the shear strains 
closely follow the yawing oscillation of the airplane. Conpeirlson of 
the load time history with the normal -acceleration record gives the 
impression that there is little or no correlation between these tvra 
quantities. The almost periodic and lightly damped characteristic of 
the loads on the vertical -tail siarface indicated by the sample records 
is also noted in the' power spectra of the vertical-tail loads for the 
whole test run shown in figirre 5* The- large power peak at a frequency 
of 0.36 cps is indicative of the response of lightly damped systems. 


Effect of Center-of -Gravity Position 

The distribution of measured peak loads on the vertical -tail surface 
for the two center-of -gravity positions shown in figure 3 cannot be 
directly compared in order to obtain the effect of center-of -gravity 
movement because the data are imcorrected for changes in turbulence 
Intensity between test runs. Inteipretlng the data for the effects of 
the center-of -gravity position is further complicated by a difference 
in moment of inertia for the two tests due to the fuel load. The gust- 
velocity spectra shown in figure 6 indicate that the turbulence inten- 
sity for the tests at the forvrard center-of -gravity position was greater 
than the turbulence intensity for the rearward center-of -gravity position 
by a factor of about 1.13 . The measured root-mean-square load for the 
forward center-of -gravity position is greater than that for the rearward 
center-of -gravity position by a factor of I. 08 . Thus, a 5-percent 
Increase in loais with rearward movement of the center of gravity and 
the simultaneous decrease in moment of inertia is indicated. 

If the calculated results shown in figure 5 were adjusted for the 
differences in turbulence intensity, the results would yield roughly 
the same root-mean-square load for both test conditions. Check calcu- 
lations showed that the decrease in moment of inertia would decrease 
the calculated root-mean-square load only about 1 percent and that the 
rearward movement of “the center of gravity (shorter tall length) would 
Increase the root-mean-square load only about 2 percent. 
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Taken Individual 1 y, ttie effect on the calculations of either the 
center-of -gravity movement or the moment-of -inertia change was veiy 
small, and the combined effects of hoth -vrould not he expected to cause 
a 5-percent difference in the root-mean-square values. In view of the 
accuracy of hoth the test results and the calculations, the indicated 
discrepancy of about 5 percent between test results and calculations is 
not too surprising. In any case, the effects on the loads in the pre- 
sent tests of the center-of -gravity movement and moment-of -inertia 
change appear very small. 


Comparison of Calculated and Measured Loads 

Figure 3 shows that the peak-load distributions on the vertical tall, 
calculated by the discrete-gust method, are considerably lower than the 
measured distributions. . The load that would be exceeded on the average 
once in 35 flight 3 niles.was taken to be 3^600 pounds. The calculated 
peak loads by the discrete-gust method for this same flight distance are 
2,150 potuads for the forward center-of -gravity position or O .60 of the 
measured load and 2,550 pounds for the rearward center-of -gravity posi- 
tion or 0.71 of the measured load. It is also to be noted that the sh^es 
of the calculated peak-load distributions based on discrete gust are con- 
siderably different from the measured distributions. For these test con- 
ditions, apparently the discrete-gust method of calculating gust loads 
seriously underestimates the gust loads on the vertical -tail surface. 

The peak-load distributions calculated from the power spectra by 
the method of reference 12 , and shown in figure 3 ^ are in better agree- 
ment with the measured results than the discrete-gust calculations, both 
in amplitude and in shape. The loads that would be exceeded on the aver- 
age once in 35 miles of flight are calculated to be 3^500 pounds for the 
forward center-of -gravity position and 2,700 pounds for the rearward 
center-of -gravity position. These values are, respectively, 0.92 and 0.75 
of the measured loads. The calculations for the forward center-of- 
gravlty positions are in good agreement with the measured values, whereas 
the results for the rearward center-of -gravity position are not nearly 
so good. 

The difference between the spectral-calculation results and meas- 
ured distributions is primarily due to the fact that the values of the 
root-mean-square load obtained from the calculated spectra were lower 
than the measured root-mean-square-load values. For the forward center 
of gravity where the calculated root-mean-square load was 9 I percent of 
the measured root-mean-square losid, good agreement was obtained between 
the calculated and the measured distributions. For the rearward center- 
of -gravity position, however, where the calculated root-mean-square load 
was only 85 percent of the measured root-mean-square load, the agreement 
between the calculated and the measured peak-load distributions was not 
so good. The 15-percent difference in root-mean-square load for the cal- 
culated and the measured spectra for the rearward center-of -gravity posi- 
tion appears to be too large for satisfactorily estimating the peak-load 
distributions . 
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It is evident ttiat estimates of the distributions of peak load, 
based on power spectra are better than the discrete-gust , estimates. For 
the present analysis, estimates of the peaks based on the spectral tech- 
niques still underestimate the loads by as much as 25 percent. This fact 
is not too surprising, however, in view of the uncertainties in the gust 
input and in the estimates of the stability derivatives. 


CONCLUSIONS 


Results of flight tests of a jet-bomber airplane flown in turbulent 
air at center-of -gravity positions of 24.3 percent and 30.6 percent mean 
aerodynamic chord indicated that turbulence excited lightly damped lat- 
ereil oscillations of the airplane. !Hie loads on the vertical tail 
resulting from the lateral oscillations were considerably greater than 
the loads predicted by discrete-gust calculations. The loads predicted 
by discrete-gust calculations were 0.60 and 0.71 of the loads resulting 
from the lateral oscillations for the test runs at forward and rearward 
center-of -gravity positions, respectively. 

Calculation of the load power spectra indicate the oscillatory 
natiire and low damping of the loads on the vertical-tall surface. The 
root-mean-square loads from the spectral calculation were in good agree- 
ment with the measTored values, being O. 9 I and O .85 of the measured root- 
mean-square load. The underestimation of the values of the root-mean- 
square load by the spectral calculations appears to be within the esti- 
mated error for the airplane stability derivatives and gust spectra. 

Calculation of the peak-load distributions from the power spectra 
by the method used in MCA Technical Note 35^0 was in fair agreement 
with the measured peak-load distributions. The eirors in the calculated 
peak-load distributions appeared to be primarily due to the underestima- 
tion of the values of the‘root-mean-sq\iare load by the calculated power 
spectra. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 8, 1956. 
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AIBPLAHa CHARACTSBIBIICB 


Oanaralt 

OrooB wine aroa, aq. ft 

Mean aen>lyiiando Ohara., ft 

Aapoot ratio . 

Llft-ourva slope for vlng, per roilan (oomiutad. froa 6A/(Avfi) vhere A is aepoot ratio) 

Wing span, ft 

Llft-ourva slope for Tertical tall per radlaa (ref. 9) 


1,173 
It. 01 


6.75 

4.63 

89.04 

1.85 


yertioal-tall gsoiiatryi 

Total Brea of vertical tall, sq ft . ■ 

7ertlcal-tall area extenaing shove horlioiital tail, eq ft 

yertioal-tall area abore etrain-gage location, sq ft . . 

Moan aarodynflalo chord, ft 

Looation of strain-gage briagai above borliontal tall, in. ....... 

Rudder and tib area, aq ft 

Aspect ratio . 

^an irtendlne ibora borlrootal tall, in 

Chord at attadhsont to borixontal tail, In 

Chord at tip, in. ' 

Fin offset, dog 

Welgjrt of vertloal-tall •truoture extending above strain-gaga bridges, Ib 


127.8 

109.4 

96.0 

9.37 

15.6 

28 

1.43 

IS 


64 


0 

360 


Hearvard oenter-of-graTlty poBltloni 

Centar-of -gravity poiltlon, percent man aerodynsirrl 0 chord 

Airplane vei^t, Ib 

Average equivalent sdrapead^ ft/aac 

Airplane naas ratio, 2H/peScg 

Oust factor, 1^ (ref. 3) 

Air denalty, p, eiuga/ou ft 

O 

Uoaant of inertia about Z-acds, slug-ft 

Tall length, ft (frtn. center of gravity to quarter chord of vortical tall) 


30.6 

60,200 

662 

22.0 

0.708 


/WOO%H 

Ve\iW.Ktr.y^ 


699,000 

52.1 


Forward oentar-of-gravlty poaltloo: 

C3anter-of-grevlty position, percent mean aarodjniunlo chord 24.5 

Airplane weight, lb 

Average equivalent airspeed, ft/soc 632 

Airplane mass ratio, 2n/paScg EL. 5 

Oust factor, Kg (ref. 3). O-TOS 

Air density, p, slugs/on ft ....... 0.002J37 

Hcoent of inertia about Z-azla, alug-ft^ 735,000 

Tall length, ft (fron center of gravity to quarter chord of vertioal tall) 32.9 
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OJABIE n 


IBEQDENCY DISOHIBtJTION OF PEAK AERODXHAMIC LOAD ON 
THE VEEITICAL TAIL IN EODGH AIR 
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TABLE IV 


STABLLITy DERIVATIVES 



(for 96 sq ft only) . 
(for 96 sq ft only) . 


Center of gravity 
Rearward Eorv/ard 

-0.25 -0.25 

0.072 0.074 

0.l4 0.15 

-0.052 -0.055 

-0.15 -0.15 

0.11 0.11 





Ruddtr. 

sngle, 

dej 




Ruddar a/7^/e 



2 |- 


v> 



A/rfi/a/fe ar/7ff/*, mecj^ured at ca/rAtr of y/\rk'/}‘y 



~rime, sec 


ro 

VO 


Figure 2 .- Ban5)le time histories. Airspeed, 390 knots} wel^t of airplane, 
60,200 pounds} zrearward center-of -gravity position. 
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Figure !<■. 
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Figure 6.- Power spectra of airspeed fluctuations. 







